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Abstract. This paper studies the queue length process in series Jackson networks with 
external input to the first station only. We show that its Markov transition probabilities 
can be written as a finite sum of non-crossing probabilities, so that questions on time- 
dependent queueing behavior are translated to questions on non-crossing probabilities. 
This makes previous work on non-crossing probabilities relevant to queueing systems and 
allows new queueing results to be established. To illustrate the latter, we prove that the 
relaxation time (i.e., the reciprocal of the 'spectral gap') of a positive recurrent system 
equals the relaxation time of an M/M/l queue with the same arrival and service rates as 
the network's bottleneck station. This resolves a conjecture of Blanc [6], which he proved 
for two queues in series. 



1. Introduction 

The queue length process {Q(t) : t > 0} of an M/M/l queue with Q(0) = admits the 
sample-path representation 

Q(t)= sup [N (t) - N^t) - (N (s) - N^s))} , 

0<s<t 

where both Nq and N\ are Poisson counting processes, iVo being the arrival process and 
N\ the process counting the number of departures and unused services. Using time-reversal 
we find that Q(t) given Q(0) = equals sup < s < t [iVo(s) — iVi(s)] in distribution, and in 
particular that ¥(Q(t) = 0|Q(0) = 0) can be expressed as a non-crossing probability 

(1) P(Q(t) = 0|Q(0) = 0) = P(T > t), 

where T = M{t > : Ni(t) < N (t)}. 

This paper generalizes (1) to a network setting, partially relying on the combinatorial 
techniques and symmetric functions discussed in [10]. The networks we study are series 
Jackson networks with external input to the first station only. We express Markov transition 
probabilities arising from these networks as a finite (weighted, signed) sum of non-crossing 
probabilities for multidimensional Poisson processes. 

Time-dependent results for queue lengths in Jackson networks are mostly focused on the 
case of two stations, e.g., Blanc [6] and Baccelli and Massey [4]. The spectral techniques from 
the work of Kroese, Scheinhardt, and Taylor [14] seem to be relevant in a time-dependent 
context, even though they are developed for the stationary distribution. O'Connell [17] 
obtains results for series Jackson networks with an arbitrary number of stations, and these 
results are most closely related to this paper. Indeed, as in [17], the combinatorial mechanism 
underlying our results is the Robinson-Schensted-Knuth (RSK) correspondence. 
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It is well-known that, through the RSK correspondence, departures from queues in series 
are related to non-crossing Poisson processes. For instance, O'Connell and Yor [18] prove 
that the cumulative departure process of the last queue is equal in distribution to the smallest 
component of Poisson processes conditioned not to cross. In turn, the latter process is 
closely related to the Charlier random-matrix ensemble. Our representations for transition 
probabilities as a sum of non-crossing probabilities have a different flavor. In fact, our 
arguments do not exploit the RSK correspondence explicitly, but we use a recent result of 
Rakos and Schiitz [20] as a starting point. Their formula is a generalization of the Schiitz 
formula for the transition probabilities in the totally asymmetric simple exclusion process 
(TASEP). 

The connection between queues in series and non-crossing Poisson processes studied here 
can possibly be used to obtain new results on queueing networks that are beyond the scope of 
the present paper. For instance, the recent work of Puchala and Rolski [19] on non-crossing 
probabilities suggests that the 'sharp' asymptotic behavior (as t — > oo) of the queue-length 
transition probabilities may be found; we discuss this in somewhat more detail at the end of 
Section 4. It may also be possible to find a 'spectral' representation for queueing probabilities 
analogous to the M/M/l spectral (or integral) representation, see, e.g., [1], from a 'spectral' 
representation for the transition probabilities of non-crossing Poisson processes. 

Convergence to stationarity. Representing Markov transition probabilities as a sum of 
non-crossing probabilities allows us to examine the speed at which the multi-dimensional 
queue length vector Q(t) converges to its steady-state vector Q(po). More precisely, given 
that all queues are initially empty, we investigate the probability that the system is empty at 
time t for large t. Using a large-deviation result for non-crossing probabilities, we determine 
the relaxation time for the queueing network, i.e., the 'asymptotic' time required for this 
probability to decrease to 1/e. For many Markovian transition probabilities (including 
reversible kernels on a finite state space), the relaxation time is the reciprocal of the spectral 
gap. We show that the network's relaxation time equals the relaxation time of an M/M/l 
queue with the same arrival and service rates as the bottleneck station, see (10). This 
resolves a conjecture of Blanc [6], which he proved for two queues in series. For work on 
relaxation times for non-Markovian processes, we refer to the recent work of Glynn, Mandjes, 
and Norros [11]. 

The model. We consider a queueing system with several single-server queues in series, 
each with an unlimited waiting space. Customers arrive according to a Poisson process at 
the first queue. They join station i + 1 after receiving service at station i, and leave the 
system after being served at the last station. The service discipline in each of the stations 
is fir st-in- first-out. All service times are exponentially distributed with a parameter only 
depending on the station, as well as mutually independent and independent of the arrival 
process. It is convenient to imagine that the arrival stream arises as the output process of 
an auxiliary zeroth station with an infinite number of customers; the service requirements at 
this station then become the interarrival times. We suppose that the total number of stations 
(excluding the auxiliary station) is N, and we let i/j be the parameter of the exponential 
service distribution at station i. In particular, vq is the arrival rate of customers at the first 
station. 

Two natural Markov processes associated to series Jackson networks are the cumulative 
number of departures vector D and the queue length vector Q. For i = 0, . . . , N and t > 0, 
let Di{t) be the total number of departures from the i-th queue in the time interval [0, t\. For 
i = 1, . . . , N, write Qi(t) for the queue length at station i at time t. The vector Q is readily 
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described in terms of D: Q(t) = (Qi(0) + Z> (*) - D^t), Q 2 (0) + D 1 (t) - D 2 (t), . . . , Q N (0) + 
D N ^(t) - D N (t)). 



Organization and notation. This paper is organized as follows. Section 2 presents results 
on the cumulative departure process which are principal tools in our investigations. Section 3 
derives two identities for the probability of an empty system at time t given an initially empty 
system, the second identity being only applicable in the positive recurrent case. Relaxation 
times are studied in Section 4. In Section 5, we discuss how the identities of Section 3 change 
when the underlying assumptions are relaxed. We review the concepts of combinatorics and 
symmetric functions that are most relevant to us in Appendix A. 

Throughout, we let the family {w n : n 6 Z) of functions on R be given by w n (t) = 
t™/n!l{ t>0 } for n > 0, while wo(t) = l{t>o} an d w n = for n < 0. Here and throughout, 1^ 
denotes the indicator function of the set A. We also define 



Unless otherwise indicated, all determinants in the present paper are (N + 1) x (N + 1), 
with the indices i and j ranging from to N. 



This section focuses on the probabilistic behavior of the cumulative number of departures 
from each of the queues, by presenting the transition kernel of this process. We also express 
this kernel in terms of the kernel arising from non-colliding Poisson processes. 

A full description of the process D on a probabilistic level is given by the transition 
probabilities <j%, defined for d,d' £ W N and t > as 



Rakos and Schiitz [20] derive an explicit formula for which we describe next. They 
study an interacting particle process known as totally asymmetric simple exclusion process 
(TASEP), which is equivalent to the departure process from queues in series. Both the 
vector-valued process D and the TASEP particle system are Markov jump processes. The 
departure process at the i-th queue Di makes a jump of +1 at rate corresponding to a 
service (or arrival if i = 0), unless Di{t) = in which case it cannot jump because the 

i-th queue is empty. Now consider the 'shifted' vector D where Di{t) = Di(t) — i, and note 
that this process has the same dynamics as the ordered positions of particles in an (N + 1)- 
particle TASEP exclusion process with particle i having jump rate vi for % = 0, . . . , N. 

The next proposition presents the result by Rakos and Schiitz in a form which closely 
follows [10]. Rakos and Schiitz [20] prove this proposition by verifying the Kolmogorov 
forward equations for D, and an alternative proof can be given using the results from [10]. 
For n£Z,t6l and < i, j < N, we write 



where (v) and h!£ 3 (u) are coefficients which are elementary functions of the vector 
v = (fo, . . . , i/jv); see Appendix A for explicit formulas. In the remainder, for notational 
convenience, we suppress the arguments v of the coefficients e£ 3 \v) and h!^ 3 \v). As an 
aside, we note that Wn 3 \t) is a shifted multiple orthogonal polynomial of Charlier [3] for 
j < i- 



W N = {x € 1 N+1 :x >...> x N }. 



2. The cumulative number of departures 



<t>t{d,d!) = F(D(t) = d'\D(Q) = d). 
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Proposition 1 (Rakos-Schiitz [20]). For d,d' G W , we have for any t > 0, 

N 

d>) = [J [e'^4~ dk ] det {w™ dj _ i+j (t)} . 

k=0 

It is our next aim to derive a representation for (f>" in terms of the transition kernel of 
Poisson processes killed on their first crossing time. To describe these processes in more de- 
tail, consider a measurable space on which a stochastic process {X(t) = (Xo(t), . . . , -Xjv(t)) : 
t > 0} is defined. Suppose that this space is equipped with a family of measures {P^ : v £ 
R^ +1 ,x G W^}. Under P£, X t is a Poisson counting process starting at Xi with rate i/j, 
and the elements of X are mutually independent. We set P^ = Pg and define T as the first 
exit time from the Weyl chamber W N , i.e., 

T = inf{/ > : X(t) W^}. 

The process X killed at time T is a Markov process on W N , and by the Karlin-McGregor 
formula [13] its transition kernel P" is given by, for t > 0, 



(3) 



A 

pnv')=p 2 ™=^T>t)=nK ; 



det < u; 



fc=0 



(*)}■ 



Note that the rates v; L do not appear in the determinant, so that a simple change-of-measure 
allows us to change jump rates: for z, z' 6 W^, 



A' 



(4) 



P-(z,z')=P t A (^^)II 



fc=0 



z k -z k 



,-(v k -\ k )t 



For z,d £ W N , we define 



A"(*,d) = ^ , " Z °---4 iV " 2iV det{^_ l 
U v (d,z) = v z Q °~ d °---v 



ZN—d,N 

N 



det 



{(-i) 



di-Zj-i+jJiN) 



"di-Zj-i+j 



As detailed in Lemma 3 below, in the terminology of Appendix A, A u (z,d) is a weighted 
analog of the number of Gelfand-Tsetlin patterns with specified bottom row and left edge. 
The representation of 4> v t in the following proposition is readily obtained from the Cauchy- 
Binet identity (e.g., [12, Prop. 2.10]) 



(5) 



^2 det {£i( z j ~ j)} det{ipj(zi - i)} = det I ^ £i(z)ipj{z) \ 
zew N kzez J 



For two kernels A and B on W N , we define the product kernel AB through AB(z,z') 

T. ye w N A i z ,y) B iy, z> )- 

Proposition 2. We have <ff = WP%k v . 
Proof. With Cauchy-Binet and the observation 



E(-i) r 4 ijv) e? = 



h 



r=0 



(-l) n eX 



(ij) 



if j < h 
if i < j, 



we obtain the claim. 



□ 
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The kernel IP is the inverse of the kernel A u as shown in [10]. Proposition 2 is therefore 
equivalent with the so-called intertwining P" A u = A u 4>^, which has the following probabilis- 
tic significance. There exist processes with values in Gelfand-Tsetlin patterns which give rise 
to couplings of the process D and the process X conditioned to stay in the Weyl chamber 
W N . The kernel A u describes the conditional distribution of D(t) given X(t) under this 
coupling. Further details can be found in [10, 22]. 

The preceding proposition connects departures from queues with non-colliding processes, 
and forms the starting point for our investigations. Another ingredient is the following 
alternative expression for A", which is Proposition 2 of [10]. It rephrases A u as a sum over 
Gelfand-Tsetlin patterns of order N + 1; see Appendix A for notation and terminology. 

Lemma 3. For z,d € W N , we have 

xeK^ : sh(x)=z, lcdgc(x)=(2 

This lemma shows that A u (z, d) vanishes unless zn = c^v, which is not so apparent from 
the definition. 

3. Transition probabilities for Q and sums of non-crossing probabilities 

This section presents two identities relating fc^(0, 0) and non-crossing probabilities, where 
k\ is the transition kernel for Q: for q, q' £ Z N , 

kUq,q')=nQ(t)=q'\Q(0)=q), 

suppressing the dependence on v on the right-hand side. 

Recall that P U (T > t) is the probability that a multidimensional Poisson process with 
rate vector v stays in the Weyl chamber W N up to time t. We write Sn+i for the symmetric 
group on {0, . . . , N}, i.e., all permutations of this set. This group acts on B^ +1 by permuting 
the coordinates. An immediate consequence of the proof of the theorem is that #(0,0) is 
symmetric in vi, . . . this has already been observed in [17]. This property does not 
extend to fcf (q, q') for q, q' / 0. 

Theorem 4. Suppose that v\, . . . , i/jv are distinct. Then we have for t > 0, 

Proof. Since = for k < and i < N, it readily follows that U u (0,z) = unless 

z = 0, in which case it equals 1. Since h^ N ^ = for k < and j < N, we similarly find 
that A u (z,d) = unless zn = d^, which can alternatively immediately be deduced from 
Lemma 3. On combining Proposition 2 with these facts, we find that for t > 0, 

#(0,0) = £>t (0,(4 ■■■>*))= Pt u (0,z)A^z, (£,...,£)) 

iez eez,zew N 



(6) = Yl P?(0,z)A. v (z,(z N ,...,z N )) 



z<=W N 



= e -£ fe =o^ £ det{w Zi - i+j (t)} Yl 

zeW N xeK^ :sh(x)=2,ledge(x)=(2 JV ,...,2 iv ) 

where the last equality follows from (3) and Lemma 3. Recall that the relevant notation is 
introduced in Appendix A. Set A v = rii<i<j<Ar(^ — v j)i a determinant of a Vandermonde 
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matrix. An empty product should be interpreted as 1. We next write the sum over x G K. N 
in terms of a Schur polynomial: 

xGK JV :sh(x)=z,ledge(x)=(zjv,...,z JV ) xe'K N ~ 1 :sh(x)=(zo-z N ,...,z N - 1 -z N ) 

= V V ■ ■ ■ V Z N S (zo-zn,-,z N -i-z n )(.V1, ■ ■ ■ , VN) 



^■■■^detj^-™^}.. 

Aj, L J jj=1,...,JV 



where the last equality follows from the alternative definition of Schur polynomials given in 
Appendix A. We thus obtain for t > 

^(0,0)= * -£"=o^ £ det^-^-W^det/if- 1 ^}.. 

The symmetry of &^(0, 0) with respect to i/i, . . . , i/jv immediately follows from the symmetry 
of the Schur polynomial. With the Karlin-McGregor formula (3) and the Leibniz formula 
applied to the second determinant, we obtain 

W, 0) = £ sgn( CT )^ 0) 1 ^ 1) 2 • • • ^ (Ar _ 2) P^)(T > t). 

" o-e5jv + i:o-(W)=0 



We next note that for a with cr(N) = 0, 

,N ,.N-1 



(-1)^ 



^— sgn(<7)i/ <y(0) i/ ff(1) ■■■^(iv_i) 

rio<j<j<Af [^(i)] _ llo<j<j<Af [^(j)] 



(-l) N sgn(a)A u n <i<j<7V-i[^(i) ~ ^(j)] 

]lo<i<j<V [ v <y{i)] Uj^WcrU) - u a{N)] _ 11^=0 ~ u <r(N)] 



Uo<i<j<N t 1 - v <?{j)/ v <?(i)\ rio<i<j<7V [ l ~ u o(j)/Vo(i)} ' 

and we have finished the proof. □ 

For j = 1, ... ,N, define the load pj of station j by pj = Vo/isj. We next work under 
the 'stability' assumption maxj^ 1 pj < 1, and express the difference k%(0, 0) — ^(0, 0) as a 
weighted sum of non-crossing probabilities. 

Theorem 5. Suppose that maxj^ pj < 1 and that all vi are distinct. Then we have for 
t>0, 

w,o)=n[i- Pi ]- e tt — nj V-" P3 ' ] / i p<TH ( r >*)- 

Proof. For A G and x e If^, we define 

u x {x)=\^---\- N XN det{\f-^}, 
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which can be expressed in terms of Schur polynomials (see Appendix A) as 



<•>*(*)= II (l-^W°---A-^(A). 

0<i<j<N A * / 

For convenience, we write oj\ for tu\(0) = rio<i<j<Af [1 ~~ ^j/^i]- 
Theorem 4 yields for t > 0, 

N 

(7) n[ 1 -fc-r 1 w,o)= E <,) p<tM ( t >*)- 

i=l o-6S JV+ i:o-(JV)=0 

Let a be the permutation for which f^(AT) < ^(jv-i) < • • • < ^a-(o)- The choice of a entails 
that the term for a converges to the constant P°"W(T = oo); the idea of the proof is to 
rewrite this term. 

As outlined in Biane et al. [5, Sec. 5.2] for the Brownian case, we have Yl {u) {T = co) = 
u)0(y)(x)\ see also [16, Sec. 4]. Moreover, u^^(-) is harmonic for the killed transition kernel 
of X under P f M, meaning that 

E H {v) (X(t) =y,T> t)u- a{v) {y) = u~ a{v) {x). 

y&W N 

Indeed, using the Karlin-McGregor formula (3), this is readily verified with the Cauchy- 
Binet identity (5). We thus obtain 

P*M(T>t) = E P Hu) (X(t)=x,T>t)[co Hl/) (x) + (l-u; a{u) (x))] 
xeW N 

(8) = W?M + E V* {v) (X(t)=x,T>t)(l-<J Hv) (x)). 

xeW N 

Next note that the expansion of the determinant uj v {x) in conjunction with the change-of- 
measure formula (4) yield 

E P"{X(t) = x,T > t)(l - u„(x)) 



N 



- - e pl w) = *,T>t) e • • • »N XN n ^« +CT(i) 

aJGVK^ aeS N+1 :cr^id 1=0 

JV 

= - E sgn(ay o ...^ oy ..^ { N N) P^(T>t) 

= - E ^; ( l)P ffM (r>t). 

<T€SjV+i:o"^id 

Replacing v by <t(z^), we see that the second term in (8) can be rewritten as 

E P^\X(t) =x,T> t)(l - u d{u) {x)) = -u a(u) E <,) pCTH (r > «)■ 
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We have thus shown that the right-hand side of (7) equals 

E ^ a(u) (T>t) 

aeS N+1 :a(N)=0 

E w ; ( i ) p^")(r>t) + i- £ ^p^)(T>t) 

E ^)P CTM (T>t) + i- E <„)P CTW (T>t) 

o-e5jv +1 :o-(Af)=0 creSjv+i 

= i- E -; ( l) pCTH (T>t), 

(765jv + i:<t(JV)^0 

and the claim follows. □ 

4. Relaxation times 

The representation in Theorem 5 is particularly suitable for studying the regime t — > oo. 
Before giving a result into this direction, we prove a large-deviation result that we need in 
our analysis. For terminology, see Dembo and Zeitouni [9]. Recall that the process X is 
defined on a measurable space equipped with P u . The P^-law of X(t)/t satisfies a large- 
deviation principle (LDP) in IR iV+1 on scale 1/t. The (good) rate function is l u : IR iV+1 — ► R 
given by 

EfcLo i x k log(x fc /i/ fc ) - x fc + "*] if x e 
oo otherwise. 

This can be seen, for instance, as a corollary of Exercise 5.2.12 in [9]. It turns out that the 
P^-law of X(t)/t also satisfies an LDP in 

Wg = {x G R N+1 : x > ... > x N } 

on the event that X is killed at the boundary of W N . 

Lemma 6. For any v G R+ +1 , the (defective) measures P u (X(t)/t G <ix, T > t) satisfy an 
LDP in on scale 1/t. The rate function is I u restricted to . 

Proof. Since P"(X(t)/t G F,T > t) < P u (X(t)/t G F), the upper bound follows from the 
LDP for the P"-law of X(t)/t, i.e., without the killing mechanism. For the lower bound, we 
prove that for any y in the interior of W$ and any 5 > 0, 

limmfilogP" (l(#e[Ife-^^*)^>*) > "^(y)- 

V k=0 J 

This inequality is trivial if y so we assume the contrary. In view of (3), the 

probability can be written as 

N 



Iu(x) 



Uk=o{ykt-st,y k t+st) J]~ LV///> 



P y (X(t) edz,T> t) 



N 

* n 

k=0 



(-) 

\VkJ 



y k t 



e -«| log(yfe/^fe)|t e — ^fei+yfei 



V 



py(X(t)/te l[( yk -e,y k + e),T>t , 



fc=0 



for any < e < 5. Since the probability on the right-hand side is bounded from below by a 
positive constant as t — ► oo, we have proven the lower bound after letting t — ► oo and then 
e -► 0. □ 
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We use the above lemma to investigate the asymptotic behavior of the sum given in 
Theorem 5, as it implies that 



(9) 



lim -logP CTM (T > t) 

t— >oo t 



By convexity, the infimum is attained on one of the faces of the polyhedron W$ . The 
following result follows from a comparison of the exponential rates of each of the terms in 
Theorem 5. We write for the i-th smallest rate among v\, . . . , z^v, so that is the 
smallest rate and i// N \ the largest. The notation f(t) ~ g(t) is shorthand for f(t)/g(t) — > 1. 



Corollary 7. Suppose that max^- =1 pj < 1 and that all are distinct. As t 

N 



oo, we have 



k?(0,0)-H[l- Pj ] 



P(i) 



n 



l<i<j<N 



[1 - v {i) / V{jl 



.p("(iV), •••,"(2),"0,''(l))£y > 



Proof. We may assume that i/q < ^i < • • • < U N since k%(0, 0) is symmetric in . . . , u^. 
Theorem 5 shows that k t (0, 0)-]jf =1 [1 — P j] is a linear combination of probabilities P* 7 ^) (T > 
i) with o-(N) / 0, and it is our aim to show that the term p^.-,^,^,^)^ > t) has the 
slowest exponential decay in t. In view of (9), we thus need to prove that a := (N, ... ,2,0,1) 
minimizes ini xeW N I a ^(x) over a 6 5at + i with <r(iV) / 0. 

Let a with a(N) / be given, and write j = <t _1 (0) < N. The rate function I a i v \ is 
strictly convex on (0, oo) N+1 . Moreover, it achieves its minimum at x = a(v) Wg\ 
since o~(v)j < o~(v)n, the latter point is separated from W$ by the hyperplane {x € M. N+1 : 
xj = xn}. The infimum of I a (v) over W$ is thus at least as large as the infimum over 
this hyperplane, with equality if and only if the minimizing argument lies in W$. It is not 
hard to see that the infimum over the hyperplane is vq + v a {v) N ~ ^^ v Q u <t{v) n i an d that it 
is achieved at v*{a) given by 



N 



otherwise. 



We conclude that, if j = o- _1 (°) < N , 



inf I a (y)(x) 
xew" 



> 



> 



inf I a{y ) (x) 

x 6 ( , oo ) N + 1 : x j = x jv 



^0 + u a(y) N 



where the last equality uses that x i— > uq + x — 2^/xvq is increasing in x > uq. 

We next argue that this lower bound is attained only for a = a. By strict convexity 
of / CT (V), the first inequality in the preceding display is strict unless v*{cr) £ , showing 
that a = a is the only candidate to attain the lower bound. In fact, since both inequalities 
become equalities if a = a, a corresponds to the unique asymptotically dominant term in 
Theorem 5. 

A direct calculation shows that for a = a, 



nf =1 [i 



-pi 



and this determines the prefactor. 



n 



l<i<j<N 



[1 - Vi/Vj] ' 



□ 
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The relaxation time is defined as the reciprocal of the exponential rate at which \k% (0, 0) — 
^oo(0> 0)| converges to zero. The proof of Corollary 7 shows in particular that, if the Vi are 
all distinct and the system is stable, the relaxation time is 

(!0) " = ; " rrr . 

^0 + ^(1) - 2 A /^p 7 (l7 v {l)K l ~ y/W)) 

This is the same as the relaxation time of an M /M / 1 queue with arrival rate vq and service 
rate v^y Therefore, the relaxation time is determined by a single bottleneck station. These 
observations prove a conjecture of Blanc [6]. 

Corollary 7 is also of interest in order to obtain exact asymptotics (as t — ► oo) for the ab- 
solute difference between fc^(0, 0) and fc^(0, 0), in contrast with the rough (i.e., logarithmic) 
asymptotics obtained here. In fact, exact tail asymptotics for non-crossing probabilities are 
a current research topic. For recent work on such asymptotics in a Brownian setting, refer 
to Puchala and Rolski [19]. Thus, a Poisson-process analog of the results in [19] would 
immediately yield new results for queues in series. 

5. Generalizations 

This section discusses analogs of the identity presented in Theorem 4 for k^{q,q r ) with 
q, q' 0, and how it changes if the service rates are not necessarily distinct. 

Define ir : W N — > Z+ via tt(x) = (xq — x\,x\ — X2, ■ ■ ■ xjv-i — xn) and recall that the 
queue length vector and departures vector are related by Q(t) = Q(0) + ir(D(t)). For £ E Z, 
we also define the 'inverse' tt£ : Z^ — ► W N of it through ne(q) = (£ + q\ + . . . + qN, £ + Q2 + 
. . . + qjy, ...,£ + q^, £), and we abbreviate ttq by tt. 

The transition probabilities for Q can be expressed in terms of those of D through 

KM) = E/eZ^te). W))> for <1,4 G Z +- Now if we define for z G W * \ Q G Z +> 

A» Q (z,q) = ^2A»(z,n £ (q))=A»(z,n ZN (q)) 

lei 

(11) = v» z °---v- N z " Y, 

xeK^: sh(x)=2, ledge(x)=7r Zjv (q) 

N-l 

= TT \ u n( - q)k ~ Zk ] det {v ZN h^ jN) ~ 1 

11 L fc J I 3 z i -z N -7r(q) j -i+3f i j=0 at_ 1 ' 

fc=0 ' J 

then it follows from P^A" = A"^ that P t " and k\ are intertwined: 

P^(z,g) = J^A^t^)) = ^A^(^^(9)) = A&fc?(z,<z). 

lei lei 

The kernel Aq is not invertible, however it does admit a (non-unique) left inverse TIq defined 
by W Q {q, z) = IF(Sr(g), z): for q, q' G Z£, 

ngA^(g,^) = 5Z H n^te)^)^^^')) = 1 {^(9)=^(?')} = l {i=i'}- 

lei zeW N lei 

Consequently we obtain the representation of the queue length transition probabilities, which 
also follows directly from Proposition 2: 

Proposition 8. For any t > 0, we have 

k\ = n%p-A%. 
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The above proposition may equivalently be written as 

(12) *4W) = E U Q(l^)K [^ Q (X(t),q'),T> t] . 

zeW N 

The function z i— > Rq(q, z) is supported on a finite set for fixed q, as is readily verified from 
its definition. Therefore, the sum over z G W N in (12) is in fact a finite sum. 

Of special interest is the case q' = 0. From (11) it then follows that Ag(z, q') can be 
expressed in terms of a Schur polynomial for nonnegative z G W N : 

xeK w : sh(x)=z, ledge(x) = (zjvr">2iv) 

It is only if the rates v , are distinct that this can be expressed as a sum of geometric terms 
and a result of the form of Theorem 4 can be obtained from (12) by suitable changes of 
measure, cf. the determinantal representation of Schur polynomials in Appendix A. At the 
other extreme, when all the Vi are equal to u, we have Aq(z, 0) = Y\o<i<j<N-i( z i ~ z i ~ 
and (12) becomes 



[] (X l (t)-A,(t)-i + i),T>t 

0<i<j<Af-l 



If g' and the rates are distinct, then Aq(z, g') has a piecewise sum-of-geometric form 
as long as z > Tr(q'), otherwise it is zero. We expect that Theorem 4 can be generalized into 
this direction, but that it involves the sum of probabilities of the form {X(t) > 7r(q'),T > t}. 

We close this section by relating Proposition 8 to the literature. In the single-station 
case AT = 1, we have A u Q (z,q) = (vi/vo) z °~ zl ~ q l{ zo > zl+q } and U^(q,z) = l{ z =( q , )} ~ 
i/ 1 l{ z= ( 9 _ 10 - ) } for q > 1 while Hq(0,z) = l{ 2= o}- This leads to 



zew 1 

z >z 1 +q' 



W zo - q (t) - ViW ZQ - q+ l(t) w zo+1 (t) 
W Zl - q -l(t) - VlW Zl - q (t) W Z1 (t) 



We know from the sum representation of the modified Bessel function (e.g., [2, 9.6.10]) that 



for dGZ. Expanding the determinant, we find that k"(q, q') equals 

p (9 '- 9)/2 ^-,+p (9 '- 9 - 1)/2 w + i + (i-p)/ E p-" 21 * 

t>q+q'+2 

where p := vq/v\ and all arguments of the Bessel functions are 2^jvQV\t. This is the well- 
known expression for the time-dependent M/M/l queue. A similar program can be followed 
for N > 1, and k\ (q, q') can then be interpreted as weighted sum of so-called lattice Bessel 
functions [4, 15, 8]. To our knowledge, for N > 2 this representation has not been recorded 
in the literature. Note that it is numerically inefficient to use the resulting representation 
of k v t , see [7]. 



12 



A. B. DIEKER AND J. WARREN 



Acknowledgments 

A.B. Dieker was financially supported in part by a postdoctoral fellowship from the IBM 
T.J. Watson Research Center, Yorktown Heights NY, USA. 

Appendix A. Symmetric functions and Gelfand-Tsetlin patterns 

This appendix defines the coefficients and used in Section 2. We also introduce 
Gelfand-Tsetlin patterns and Schur polynomials, which play important roles in the proof of 
our main results. More details can be found in, e.g., Chapter 7 of Stanley [21]. 

The rth complete homogeneous symmetric polynomials in the indeterminates ceo, • • • &n 
is given by 

h r (a)= Yl 

fco>0,...,fcjv>0:fco+fciH \-kfj=r 

By convention ho = 1 and h r = for r < 0. Now for < % < j < N, let h { j j \a) = h r {a {i ^) 
where afa') is the (N + l)-vector (0, . . . , 0, a^+i, aj+2> • • • > a ji 0, • • • 0) obtained from a by 
setting the first i + 1 weights and the last N — j weights equal to 0. Equivalently it is the 
rth complete homogeneous symmetric polynomial in the indeterminates ctj+i, . . . ,ctj. We 

set hi jj \a) = l(r = 0). 

We write e r for the rth elementary symmetric polynomial defined as 

e r (a) = X a kl ---a kr . 

0<fcl<fc 2 <---<fcr<A r 

In analogy with the complete homogeneous symmetric polynomials, we use the conventions 
e^\a) = l(r = 0) and eo(a) = 1. We also set e^\a) = e r (a^), so that in particular 
e { r j \a) = for r < 0. 

Let x be an array of real-valued variables x = (x°, . . . , x N ) with x k = (x*, x\, . . . , xfy 6 
Z fc+1 , such that the coordinates satisfy the inequalities 

~k < k-l < k < k-l < < k < k-1 < k 
•^k — x k~l — x k~l — x k-2 — • ■ ■ — x l — x — x 

for k = 1, . . . , N. We write K w for the set of all x satisfying the above constraint, and 
say that any x G K w is a Gelfand-Tsetlin (GT) pattern of order N + 1. For x G K. N , we 
set sh(x) = (xq ,Xi , . . . , x^) and ledge(x) = (xq, . . . ,x^). If ledge(x) > 0, i.e., all xj are 
nonnegative, integer-valued GT patterns of order N + 1 are in one-to-one correspondence 
with so-called semistandard Young tableau with N+ 1 rows and entries not exceeding N+ 1. 
For a vector a = (cyq, cci, • • • , ajv) of weights, we define the weight q x of a GT pattern x by 

N 

n 

k=i 

The Schur polynomial in the indeterminates ao, • • • ckjv corresponding to a nonnegative 
z G PF^ is given by 

xeK iv :sh(x)=z 

Schur polynomials are symmetric in the which is readily seen from their alternative 
determinantal definition: 



« x = V 11 a k 



s z (a) = r det {a^ l+7V ) . 

I Ii> , V "' " /\ 1 ' J 
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This representation requires that the a« be distinct, but there is no singularity; there are 
factors of (ctj — ay) implicit in the determinant which cancel with factors in the denominator. 
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